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TECHNICAL MEMORANDUM X- 64948 


THE MOTION OF THROW AWAY DETECTORS RELATIVE 
TO THE SPACE SHUTTLE 

I. INTRODUCTION 


This report documents the results of a study of the motion of TAD’s relative to 
the Space Shuttle. It was done in support of the AMPS program. A TAD is a Throw 
Away Detector — a small, inexpensive, expendable satellite that has no guidance, control, 
or navigation capability and no propulsive capability. A TAD may be spin stabilized. The 
TAD’s will probably be ejected from the Shuttle by a spring device with a very small 
impulsive velocity on the order of 1 m/s. 

The TAD’s may be utilized as diagnostic tools in the AMPS program to examine 
the electromagnetic interference environment around the Shuttle and to sample the gas 
contaminants emitted by the Shuttle. They may also be used to measure the perturbations 
induced by the Shuttle in the ionosphere, i.e., for plasma wake and sheath measurements. 
There are probably other applications also. The maximum range expected to be utilized 
by the TAD’s will probably be approximately 10 km and in many cases the region of 
interest may be less than 1 km from the Shuttle. 

Once the TAD’s are ejected from the Shuttle, the experimenter has no more con- 
trol over their motion; thus, the only control over their motion is through the ejection 
conditions. Those conditions must be chosen so as to make the TAD’s execute the 
desired motion to the extent possible. The ejection must also be designed to avoid recon- 
tact between the TAD and the Shuttle. 

In this report, we shall be concerned with determining: (1) what relative motions 
are possible with the TAD’s, (2) the sensitivity of the motion to the ejection conditions, 

(3) how to avoid recontact, and (4) how the drag force affects the relative motion. We 
shall not be concerned with hardware, e.g., with designing a spring ejection device, nor 
shall we be concerned with communications problems between the TAD’s and the Shuttle. 
These problems must be covered under separate studies. 

A method is developed that enables one to target a TAD to a specific point at a 
given time. A method is also developed that allows one to map specified allowable errors 
at the target point back into allowable errors in the ejection conditions. This allows one 
to determine how accurately the ejection mechanism must be built. 

It is shown that the motion of a TAD thrown out on the left side of the Shuttle 
is the “mirror-image” of the motion of a TAD thrown out on the right side; thus, the 
analysis of the motion has to be made only for ejections on one side. 


The targeting capability allows one to examine whatever point or points may be 
of interest. The only limitation to the number of points examined will be the number of 
TAD’s available because, in general, it will require a separate TAD for each individual 
point (unless the TAD’s are recovered and reused). 

Finally, the effects that an uncertainty in the differential drag might have on the 
targeting problem are examined. It is demonstrated that the effects are small at high 
altitudes but not at low altitudes. 

Most of the results are presented in rather general or parametric form. This was 
necessary because the AMPS program is still in a formative stage and no specific applica- 
tion has yet been defined. It is believed that the techniques presented here are, for the 
most part, developed to the extent that they may be directly applied to any specific case. 


II. FORMULATION OF THE PROBLEM: THE MOTION OF ONE FREE 
PARTICLE RELATIVE TO ANOTHER IN AN INVERSE SQUARE 
GRAVITY FIELD IN THE PRESENCE OF ATMOSPHERIC DRAG 


In this section the problem of relative motion, i.e., the motion of one orbiting 
particle relative to another, is formulated and solved by making certain approximations 
such as neglecting second-order small quantities. There is nothing new about the develop- 
ment; it has been known for at least 15 years [1]. It is included in this report in order 
to make it self-contained and to make clear the assumptions made in arriving at the 
solutions. 

Consider a Shuttle orbiting the Earth in a circular orbit of radius R. Establish a 
local coordinate system at the center of mass of the Shuttle with y positive up along 
the local vertical, x in the orbit plane with the positive direction opposite the orbital 
velocity vector, and z perpendicular to the orbit plane and positive in the direction of 
the orbital angular momentum vector as shown in Figure 1. A TAD may be located at 

/? relative to the center of the Earth and at r relative to the Shuttle. 

By definition, the vectors r , R, and p are related by 


->• -»■ 
p = R + r 


( 2 - 1 ) 


Clearly, the coordinate system just described is a rotating one, rotating at the orbital rate 
w of the Shuttle. The rotational rate co is related to the orbital radius R by the 
classical two-body result 


cu 2 R 3 = 




( 2 - 2 ) 



where p is the product of the universal gravitational constant and the mass of the Earth. 
This relation will prove useful shortly. 


To formulate the equations of motion of a particle, an expression for the accelera- 
tion of its position vector must be obtained and then, according to the Newtonian scheme, 
this must be equated to the forces acting on the particle divided by the mass of the par- 
ticle. From fundamental considerations [2] the time derivative of a vector in an inertial 
system is related to the time derivative of the same vector in a rotating system by the 
relation 


■■■ y — ^ — v 

P = [p] + w X p = 



I 

+ co X) p 


(2-3) 


where p represents the velocity of the particle in the inertial system, [p ] represents the 

velocity of the particle in the rotating system, and co is the angular velocity of the 
rotating system. In this application the orbital rate of the Shuttle is given by equation 
(2-2). A second application of the operator listed in equation (2-3) gives the acceleration 

of p ; namely, 


p" = [p ] + [co ] X p + 2 co X [p ] 4' co X (co X p) . (2-4) 


The second, third, and fourth terms on the right side of equation (2-4) are referred to as 
the transverse force, the coriolis force, and the centrifugal force, respectively. 

...y t t 

If the Shuttle is moving in a circular orbit, co is constant and co is zero. Tins 
makes the second term on the right side of equation (2-4) vanish. We shall make this 
assumption for this study. 

Equation (2-4) is expressed in terms of ~p . For the relative motion analyses it 

needs to be expressed in terms of ~r . This is easily accomplished by using equation 
(2-1). Differentiating equation (2-1) with respect to time yields 

p = t + t . (2-5) 


Applying the operator listed in equation (2-3) to all terms yields 



Since the Shuttle is in a circular orbit, R is constant so that [R] , the velocity of R in 
the rotating system, is zero. The second and fourth terms on the right side of equation 
(2-6) can be combined by application of equation (2-1), and the combined term then 
cancels the second term on the left side to finally yield 

[pl = [t] . (2-7) 


Similarly, it can be shown that 

[?]=[?] • ( 2 - 8 ) 


With these results, equation (2-4) can be rewritten as 
p = [?] + X [tj + Z X (w X p) 


(2-9) 


By Newton’s Second Law, the inertial acceleration of the TAD, p , can be replaced by 
the force per unit mass acting on the TAD. Two forces are considered, the inverse square 
Newtonian gravity and the drag due to the atmosphere. Some rather tortuous logic is 
required to arrive at the next point, which is that the force acting on the TAD is the 
inverse square gravitational force plus the differential drag force or the difference in the 
drag force between the TAD and the Shuttle. 

Differentiating equation (2-1) twice yields 

P = R + r" (2-10) 


and by Newton’s Second Law 

(F/m) TAD = (F/m) Shutt i e + A(F/m) . (2-11) 

The interpretation of the last term is that if there is an acceleration of the relative posi- 
tion between the TAD and the Shuttle, there must be a differential force between them. 
Unit masses can be considered for the moment to eliminate m from equation (2-11), 
and the forces originating from gravity and drag can be written as 


4 


( 2 - 12 ) 



Now R and co can be approximated as constants which means, in effect, that the 
effects of drag on the Shuttle are ignored; hence, the second term on the right side of 
equation (2-12) is dropped. It is not, however, ignored in the last term; instead the last 
two terms are written as 


— ^ 


F 


l D 





af d 


D 


(2-13) 


The justification for dropping Fg^ in one place and leaving it in another is that the 

net result of the drag force on one vehicle for a short time may be small but the result 
of the difference in drag between two vehicles may be significant. 

The first and fourth terms on the right side of equation (2-1 2) cancel leaving 


F t +F t =F t + AF d = F t +D . (2-14) 


This is the desired result; namely, that the force on the TAD is the inverse square gravity 
force plus the differential drag force. Thus, we can write 

jl — y — > . 

p = F t + D = ~ + D . (2-15) 

g P 


This is substituted into equation (2-9) which is then solved for [ r ] . The result is 


[t] = ^ + D-2^X [t] -w X(w Xp) 
P 3 


(2-16) 


It can be verified, by direct multiplication, that the last two terms in equation (2-16) 

have no z-component. Only the first term on the right side has a z-component (since 
— 

the differential drag D is taken to be opposite the velocity vector and, therefore, along 
x). This can be equated to the z-component on the left side to obtain 



5 



Equation (2-16) can be written again as it stands, but with no z-component in any term; 
it must now be interpreted as being two dimensional 



- + D - 2co X [t] -2X fcoXp) 

P 3 


(2-18) 


In two dimensions the last term can be written as co 2 p . This can be verified by direct 
multiplication. This can be combined with the first term, where p can be replaced by 
equation (2-2), to yield 




+ D - 2co X 



(2-19) 


An expression is now developed for 1 jp . By definition 


->■ 

P 


-> r* 
r + R 


so that 


p 2 = p ■ p = (7 + R) • (r + R) = r 2 + 2rR cos 9' + R 2 . (2-20) 


This last relation can be verified by referring to Figure 1. From equation (2-20) one can 
get 


P 



(-cos 6') + 



( 2 - 21 ) 


The radical is recognized as the generating function for the Legendre polynomials, 

P n (-cos 9'), so that 

oo 

W 2 p » ( - cosO,) (0 n • (M2) 

n-0 


6 


i 


The ratio (r/R) is a small quantity, so to first order in this small quantity. 


i = 1 1 1 . cos fl - (i)j 


(2-23) 


If the cube of equation (2-23) is taken and first-order small quantities are retained, one 
then obtains 


7 = gj [' ■ 3 cos 01 (0j 


(2-24) 


This result substituted into equation (2-19) will yield 


r = c p 


(3 cos d'.0 + D - 2co X [7] 


(2-25) 


By inspection of Figure 1 it can be seen that 


r cos 6' = y 


so that equation (2-25) becomes 


7 = 3co 2 y £ + D- 2uX [7] 

R 


(2-26) 


Since p is now two-dimensional, as it has been since equation (2-18), one can write 


7 _ lx + jty + R) 


R 


R 




(2-27) 


Again x/R and y/R are very small quantities compared to 1 (about 1CT 3 ) so one can 

approximate p^/R by t so that equation (2-26) becomes [this approximation is very 
important to the further development of the equations since it linearizes them] 


[r”] = 3oj 2 y6 + D -2w X [7] 


(2-28) 


7 



D can be taken to be of constant magnitude and opposite the velocity vector, i.e., in the 
+x direction. Multiplying out the last term yields 


ix + jy = (3co 2 y)f + (D)i + (2ccy)i + (-2ccx)f 


(2-29) 


Equating components of the unit vectors yields 


x = D + 2coy 
y = 3co 2 y - 2ccx 


(2-30) 


These are two second-order linear differential equations with constant coefficients that 
can be solved with relative ease. They are coupled in the sense that the x-motion affects 
the y-motion and vice versa. 

In the z-equation (2-17) one can replace p by co 2 R 3 from equation (2-2) and 
then approximate (R/p) 3 by 1 to get 


z + w 2 z = 0 


(2-31) 


The z-motion is simply that of a harmonic oscillator with the frequency equal to the 
orbital rate; i.e., it completes one oscillation in one orbital period. In arriving at the 
z-equation a first-order small quantity (y/R) was ignored. The same was done in equation 
(2-27). These made possible the linearization of the differential equations that are then 
solvable in closed form. 

A great deal of experience with these equations indicates that they are quite good 
at least for short time periods or where the TAD’s do not drift too far away from the 
Shuttle in the vertical direction. Most of the neglected terms are at least 10" 3 times the 
terms that were retained. 

The great value of these equations is that no numerical integration is required to 
analyze the motion of the TAD’s. In summary, the following assumptions were made to 
arrive at these equations: 

• The Shuttle is in a circular orbit (w =0). 

• The Shuttle orbit is unaffected by drag (co and R remain constant). 

• The differential drag force D is constant and opposite V . 


8 



• The gravity field is inverse square (oblateness is neglected). 

• Second-order small quantities are neglected in the expansion of l/p . 

• First-order small quantities (1 0" 3 compared to 1) are neglected in linearizing 
the differential equations. 

If any of these assumptions should be grossly violated, the resulting equations, i.e. (2-30) 
and (2-31), most likely would be invalidated. For example, if the Shuttle should be in a 

highly elliptical orbit, then co would not be zero and could not be ignored in equation 
(2-4). Other violations of the assumptions can be easily determined. 

The solutions of these equations are obtained rather easily. The solution to the 
z-equation (2-31) can be written as 


z = A cos cot + B sin cot (2-32) 

and if at t = 0 , z =z 0 and z = z 0 , the constants A and B are A = z 0 and B = 
z 0 /co ; therefore, the solution is 


z = z 0 cos cot + z 0 /co sin cot 
z = -coz 0 sin cot + z 0 cos cot 


(2-33) 


Equations (2-30) can be solved in various ways. In the method chosen here they 
are converted to four first-order equations by the following change of variables: 


Xj = x 


x 2 = y 
x 3 = x 
x 4 = y 


(2-34) 


The differential equations then become 


9 






ffl 


H 


or compactly as 


x = Ax + B 


with the solution [3] 


x(t) = e At x(t 0 ) + j e" 


A is, of course, the 4X4 matrix defined in equation (2-35) and e Al; is defined as 


e At = 1 + At + ~~~ 2 + . . 

ot 


For simple matrices, and in this case in particular, equation (2-38) can be used as a com 
putational technique for e At . By direct multiplication and summing of the elements, 
e A ^ can be verified to be 


,At = 


1 6cot - 6 sin cot 4/co sin cot - 3t 2/co(l - cos cot) 

0 4-3 cos cot -2/to(l - cos cot) 1/co sin cot 


6co(l - cos cot) 4 cos cot - 3 


3co sin cot 


-2 sin cot 


cos cot 


Then e"^ can be obtained simply by replacing t by -t . If t 0 is chosen as zero and 
all of the manipulations implied in equation (2-37) are carried out, the solutions can be 
verified to be 


3 

x = x 0 + c 0 t - - Dt 2 - 2b 0 sin cot + 2a 0 (i - cos cot) 

2c 0 2Dt , , , . , 

y = — b 0 cos cot + a 0 sin cot 

3co co 

x = c 0 - 3Dt - 2cob 0 cos cot + 2coa 0 sin cot 
2D , 

y = + cob 0 sin cot + coa 0 cos cot , 

CO 


where 


a = y° + 2D 
a 0 - — ■ t — 

co co^ 


b 0 = 3y 0 - 


2x 0 

co 


c 0 = 6coy 0 - 3x 0 


(2-40) 


(2-41) 


Equations (2-33) and (2-40) combined with the initial conditions (x 0 , y 0 , z 0) x<), y 05 z 0 , 
co, and D) completely determine the motion of a TAD relative to the Shuttle. 


III. INITIALIZATION OF THE MOTION: THE ORBITAL RATE, 
THE DIFFERENTIAL DRAG, AND THE 
EJECTION CONDITIONS 


The orbital rate co appearing in the relative motion equations depends only on 
the orbital radius R of the Shuttle, which is the radius of the earth Rg plus the altitude 
of the Shuttle orbit. It is given by equation (2-2) or 


co = \J ill R 3 


> 


(3-1) 


where m = 3.986012 X 10 14 m 3 /s 2 . 
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The differential drag D appearing in the relative motion equations is the differ- 
ence in drag between the TAD and the Shuttle and is given by 


D = ^pv 2 


c D a> 


m 


TAD 


CdA\ 

m /Shuttle 


(3-2) 


In this equation p is the atmospheric density and v is the orbital speed, assumed here 
to be the same for the Shuttle and TAD since they are in nearly the same orbit. For 
circular orbits, there exists the simple result 


v 


2 



(3-3) 


The quantity (Cj-jA/m) is the ballistic coefficient of the object, Cq is the coefficient of 

drag with a typical value of approximately 2, A is the effective cross-sectional area pre- 
sented to the resisting medium, and m is the mass of the object. From equation (3-2), 
it is seen that a differential drag exists only if there is a difference in the ballistic coeffi- 
cients of the TAD and the Shuttle. 

The speed used in equation (3-2) should be the differential speed between the 
orbiting vehicle and the rotating atmosphere. The error caused by using the inertial speed 
however is slight and can be offset if necessary by “fudging” some other parameter, for 
example Cp or A . 

The ballistic coefficient of the Shuttle depends on its attitude and on the altitude 
of its orbit. At altitudes near 400 km the ballistic coefficient varies in the range 

0.002 ^ < ffehuttle < 0-009 g . (3-4) 


because of attitude variations. For this study, the TAD’s are represented by 45.7 cm 
(18 in.) spheres weighing 22.7 kg (50 lb). The effective cross-sectional area is assumed 
to be a circle of 45.7 cm (18 in.) diameter, giving A = 0.164 m 2 . The coefficient of 
drag is taken to be 2. For the ballistic coefficient, this gives 


^TAD 



= 0.0145 — 
kg 


TAD 


(3-5) 



An intermediate value of the ballistic coefficient of the Shuttle between the two limits 
listed in equation (3-4) is taken. The intermediate value is taken as 0.0045 m 2 /kg so 
that the difference is a convenient number for computations, 


CdA\ /c d a\ 

111 /TAD \ rn / Shuttle 


0.0145 


nr 

kg 


- 0.0045 


nr 

kg 


- 0.01 


m 

kg 


(3-6) 


At 400 km the atmospheric density, according to the 1962 U. S. Standard 
Atmosphere, is 


p 400 — 6.5 X 10- 12k -S 
nr 


(3-7) 


The square of the orbital speed at 400 km altitude is 


3.986012 X 10 14 ™ 


a = 

R 


6.778160 X 10 6 m 


s2 = 5.88067 X 10 7 ~ 

s 2 


With these values, the differential drag from equation (3-2) is 


(3-8) 


D = 1.9 X 10' 6 ^ . (3-9) 

s 2 

Variations in the attitude of the Shuttle can make this value increase or decrease. Using 
the maximum and minimum values of the Shuttle ballistic coefficient gives 


1.05 X 10" 6 — <D < 2.4 X 10' 6 ^ . (3-10) 

s 2 s 2 

These values could probably be altered considerably by different designs of the TAD’s, 
i.e., by varying their size, shape, and weight. Operating at different altitudes will also 
vary D through its dependence on atmospheric density. 

If the constant difference of ballistic coefficients of the TAD and Shuttle of 
0.01 m 2 /kg is assumed, then D varies with altitude in the following way: 
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Altitude (km) 

p(kg/m 3 ) 

v 2 (m 2 /s 2 ) 

D(m/s 2 ) 

200 

3.3 X 10' 10 

6.05946 X 10 7 

1.0 X 10' 4 

300 

3.6 X 10' 11 

5.96873 X 10 7 

1.07 X 10' 5 

400 

6.5 X 10’ 12 

5.88067 X 10 7 

1.91 X 10’ 6 

500 

1.6 X 10" 12 

5.79517 X 10 7 

4.6 X 10' 7 

600 

4.6 X 10' 13 

5.71212 X 10 7 

1.3 X 10’ 7 

700 

1.5 X 10' 13 

5.63142 X 10 7 

4.2 X 10" s 


A plot of the differential drag versus altitude is shown in Figure 2. These values should 
be taken as order of magnitude values. 

It is assumed that the TAD’s are carried in the Shuttle bay until some point in 
the mission, at which time they are ejected impulsively with arbitrary speed and in an 
arbitrary direction relative to the coordinate system shown in Figure 1 . Because of a 
lack of knowledge of precisely where the TAD’s may be located in the Shuttle, it is 
assumed for this study that they are at the center of mass of the Shuttle, i.e., at the 
origin of the relative coordinate system. This defines the initial position for the relative 
motion studies as 


x 0 = 0 

y 0 = 0 

z 0 = 0 (3-11) 


The direction and magnitude of the initial velocity can be arbitrary. Consider Figure 3 
to initialize the velocities. 

If the TAD’s are ejected from the Shuttle in an arbitrary direction with an arbi- 
trary impulse, then from Figure 3 it can be seen that the initial velocities are given by 


x 0 = AV cos 0 sin tp 
y 0 = AV sin 0 

z 0 = A V cos 0 cos ip . (3-1 2) 
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These quantities (co, D, x 0 , y 0) z 0 , x 0) y 0 , and z 0 ) completely determine the relative 
motion. 


IV. ANALYSIS OF THE MOTION 


The equations of motion (2-33) and (2-40) with the initial conditions defined in 
Section III are 


x = -3(AV cos 9 sin <p)t - — Dt 2 + — (AV cos 9 sin (/?) sin cot + 2 ( AY 7 n . ^ 

2 co V co 


co 


') 


+ ^5 ) (1 - cos cot) 


CO 


CO 


2 a ■ i 2Dt , / AV sin 0 2D , . , 

y = (AV cos 0 sm <p)(l - cos cot) + I + — ) sm cot 


co co 


z = — (AV cos 0 cos i p) sin cot 

CO 


(4-1) 


For purposes of discussion a 400 km altitude circular orbit is assumed with 


co = yj ^ ~ 1-13 X 10" a (radians/s) 


and 


-i- = 884 s 
co 

The differential drag D is taken to be 
D = 1 X 10' 6 m/s 2 


The z-motion is that of a simple harmonic oscillator with a period equal to the orbital 
period. The amplitude of the motion is 

AV a 

z max = 77 CQS ° C0Sip 
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If the TAD was kicked out in the z-direction ( 9=0 degree, = 0 degree) with an 
impulse of 1 m/s, the maximum amplitude would be 884 m; if the impulse was 10 m/s, 
the amplitude would be 8.84 km; and if the impulse was 0.1 m/s, the amplitude would be 
88.4 m, etc. If the TAD was kicked striaght up or straight down ( 9 = 90 degrees or 
270 degrees) there would be no z-motion. If the TAD was ejected anywhere in the x-y 
plane (Fig. 3) where p - 90 degrees or 270 degrees there would be no z-motion. 

The y-equation, the one determining the vertical motion, has a secular term in it, 
-2Dt/co. At the altitude considered here, the value of 2D/co is approximately 2 X 10" 3 
m/s. This gives a decay rate for the TAD relative to the Shuttle of about 7 m/hour. At 
an altitude of about 300 km where the drag is an order of magnitude higher, the decay 
rate would be about 70 m/hour, and at an altitude of 200 km, the decay rate would be 
about 700 m/hour. 

The first and third terms in the y-equation are periodic terms giving rise to an 
oscillatory motion in the vertical direction. The amplitude of the first term at one-half 
period is (4/co)AV cos 9 sin p. If 0 = 0 degree, p = 90 degrees, and AV = 1 m/s, this 
term would be approximately 4 km. This amplitude can be suppressed by choosing proper 
values of 9 and i p. The amplitude of the last term is (AV/t o) sin 9 which for 9-9 0 
degrees and AV = 1 m/s is approximately 1 km. This too can be suppressed by keeping 
9 sufficiently far away from 90 degrees or 270 degrees. 

The first two terms in the x -equation are secular terms that will monotonically 
increase or decrease the x-coordinate, i.e., the coordinate in the direction of orbital 
travel. This will make the TAD get progressively further in front or progressively further 
behind the Shuttle with increasing time. For small values of t the first term completely 
dominates the second. Equating the terms enables one to find the approximate time at 
which they become of the same order of magnitude. This gives 


2AV cos 9 sin p 
D 


For AV = 1 m/s, D = 10' 6 m/s 2 , and 2 cos 9 sin < p- 1, the above equation gives about 
10 days. This means that for most applications, where we are interested in the motion 
for at most a few hours, the first term is completely dominant. The second term is 
always in the -x direction or in the direction of the velocity vector. Thus, after it begins 
to dominate the first term, the x-coordinate will become progressively more negative and 
the TAD will move out in front of the Shuttle. 

The sign of the first term depends on the ejection conditions and can be positive 
or negative. Thus for the first few days of the motion the x-coordinate can be made 
positive or negative by proper choice of the initial conditions. The magnitude of the first 
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term can vary between ±3AV which, for a separation velocity of 1 m/s, is ±3 m/s. In 
1 min (60 s), this term can be as large as ±180 m. In one orbital period (90 min or 5400 
s) this term can be as large as ±16 km. Thus, this term can lead to “run-away” motion 
in the x-direction if the ejection conditions are not carefully controlled. 

If it is desired to keep the TAD within a certain distance of the Shuttle for a 
certain length of time, for example keep x less than x max for a period of time t max , 

the following conditions must be imposed on the ejection velocity 


a \ r n ^ x max 

AV cos 0 sin <p < 

3t 

J L max 

As a numerical application, suppose that it was desired to keep x less than 1 km for one 
orbital period. The above relation would yield the maximum value for AV cos 0 sin ip 
to be 0.061 m/s. If the ejection velocity was 1 m/s, the product cos 0 sin ip would have 
to be 0.061. For i p near 90 degrees or 270 degrees (in-plane ejection), 0 must be 
near 90 degrees or 270 degrees (vertical ejection). For ip near 0 degree or 180 degrees 
(sidewise ejection), 0 could be arbitrary in value. 

In summary, to keep the x secular motion in bounds the ejection velocity must 
be approximately perpendicular to the orbital velocity direction. This is in any direction 
in or near the y-z plane. If it was required to keep cos 0 sin ip < 0.061 and if 0 = 0 
degree, the maximum value of ip would be 3.5 degrees. 

The third and fourth terms in the x-equation are oscillatory with a period equal to 
the orbital period and having the following amplitudes, respectively, 


4AV c os 0 sin ip 
co 


and 


4AV sin 0 

GO 


For separation velocities of 1 m/s these amplitudes can be as large as 4 km. The ejection 
conditions must be designed to keep them to whatever magnitude is desirable. If the 
product AV cos 0 sin i p was kept small to limit the x secular motion, the first ampli- 
tude would also be small. Using the previous example where the product AV cos 0 sin <p 
was limited to 0.061, the first amplitude listed above would be limited to less than 250 m. 
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The product AV cos 0 sin ip could be kept small by choosing \p small. However, 
6 could still be near 90 degrees and the second amplitude could be as large as 4 km (for 
separation velocity of 1 m/s). If it is likewise desirable to keep the amplitude of the 
second term below a certain value x max , the following condition is imposed on 0 


0 < sin -1 


6jx max \ 

, 4AV / 


For x max = 1 km and AV = 


m/s this gives 


0 < sin" 


^io ' 3 • io +3 j 


< sin -1 (0.25) 


or 


0 < 15 ° 


To satisfy all of these conditions means that the ejection must be nearly sidewise; 

< p within ±3.5 degrees of 0 degree or 180 degrees and not far above the local horizontal, 
0 < ±15 degrees. This condition results from a separation velocity of 1 m/s with the 
restriction of limiting the motion of the TAD to within a radius of 1 km of the Shuttle 
for a time of at least one orbital period (90 min). 

The first term in the x-equation determines whether the TAD gets progressively 
in front of or behind the Shuttle. If x is negative, the TAD is in front of the Shuttle; 
if x is positive, the TAD is behind the Shuttle. If the ejection is in the positive z- 
direction, then 0 is small, for example 0 < ±15 degrees, and cos 0 is a positive num- 
ber near 1. <p can be a small positive or negative angle, <p < ±3.5 degrees, and sin i p 
will likewise be a small number, positive if < p is positive and negative if \p is negative. 

If i p is positive (backward ejection), the whole first term is negative and if <p is nega- 
tive (forward ejection), the whole first term is positive. Thus, if it is desired to put the 
TAD behind the Spacelab, the ejection must be forward and if it is desired to put the 
TAD in front of the Spacelab the ejection must be backward. 

The coefficients in the y-equation are the same as those in the x-equation so that 
if the x-motion is kept small, the y-motion will likewise be kept small. The coefficient 
hi the z-equation, however, is different so that it is large for small 0 and small <p . The 
only way then to keep the z-motion small would be to reduce the impulse (AV) given to 
the TAD. An impulse of 1 m/s, however, keeps the z-motion to less than 1 km. If it is 
desired to keep it to approximately 100 m, the impulse would have to be reduced to 
about 0.1 m/s. 
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Another very important aspect to consider in the analysis of the motion of the 
TAD’s is collision or recontact with the Space Shuttle at some time after the ejection. 

It would seem that these should be avoided even though recontact would be at a relatively 
slow velocity, i.e., at about whatever velocity the TAD was ejected. A collision is defined 
as the condition where all of the coordinates (x, y, and z) become simultaneously zero or 
near zero at some time after ejection. Collisions are avoided by avoiding the ejection 
conditions that result in this occurrence. 

To avoid “run-away” x secular motion or large amplitude x and y oscillatory 
motion, the ejection will have to be in the sidewise direction as described earlier. For 
sidewise ejections there is always a z-motion and it has zeroes only at t = T/2, T, 3T/2, 
2T, . . . etc. Thus, to examine the collision possibilities, it is necessary to examine the x 
and y motions only at half and whole integral periods. 

At half-integral periods, cot = 7 r, 37 t, . . . etc., and 


, T ft7T , _ , 

t = n — = — - ; cos cot = -1 

2 co 


sin cot = 0 


At these times the x, y, and z coordinates from equations (4-1) are: 


x 



= -3(AV cos 0 sin i p) 


mr 

co 


3 p n 2 7 T 2 + 4 / AV sin 0 
2 co 2 \ co 




4 , AX7 n . v 2D nir 

— (AV cos 0 sin 14 ) — 

co co co 



n = 1, 2, 3 . . . 


(4-2) 


To find the collision possibilities, these expressions are simultaneously equated to zero. 
This, then, places certain requirements on the ejection conditions. At a half-integral period 
y equal to zero requires 


AV cos 0 sin = - 


mrD 

2co 


n - 1, 3, 5, . . . 


(4-3) 


For small n (for example 1, 3, or 5) and for altitudes around 400 km where D ~ 

10 ‘ 6 m/s 2 , the order of magnitude of the right side ranges from 10 ' 3 to 10‘ 2 . For AV’s 
of approximately 1 m/s, the product cos Q sin must range from 10 “ 3 to 10" 2 . Since 


19 


we are dealing with nearly sidewise ejections where 9 < 15 degrees, cos 6 is near 1, so 
that for the product to be small sin </? must be small. If sin i p ranges from 10~ 3 to 
10" 2 , ranges from 0.05 degree to 0.5 degree. Thus, ip should be kept larger than 
0.5 degree to keep y from becoming zero at half-integral periods. 

If the condition (4-3) is substituted into the x-equation of (4-2), the first two 
terms cancel exactly. If the remaining expression is equated to zero, the resulting condi- 
tion is 


AV sin 0 


2D 

co 


(4-4) 


The right side here is of the order of magnitude 10' 3 . For a AV of 1 m/s this requires 
that 0 be very near 0 degree or 180 degrees (horizontal ejection). 

Thus, to avoid collisions at half-integral periods, one must avoid exactly sidewise 
ejections, i.e., straight out the z or -z directions. The direction of the ejection should 
be at least 0.5 degree from the z-axis, i.e., one should have |ipi > 0.5 degree. 

At whole integral periods, cot = 2ir, 4 tt, . . . etc., and 


t 


nT = 


2m 

co 


cos cot = +1 


sin cot = 0 


At these times the x, y, and z coordinates from equations (4-1) are 


x(nT) = -3(AV cos 9 sin ip) 


2m 

co 


3 D 4n 2 7T 2 
2 co 2 


y(nT) = - — 


2D 2n7r 


co co 


z(nT) = 0 n = 1, 2, 3, . . . 


(4-5) 


The y-equation shows that there is no control, through the ejection conditions, 
over the y-coordinate at whole orbital periods. In fact, if there is no differential drag, 
the y-coordinate comes back exactly to zero at every orbital period. At one whole orbital 
period (n = 1) at an altitude of approximately 400 km where D ~ 10" 6 m/s 2 , y is 
approximately -10 m. This is uncomfortably close to zero. Collisions will have to be 
avoided by control of the x-motion. Equating x to zero give? the following condition . 
for the ejection 
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(4-6) 


. 117 rD 

AV cos 0 sin ip = 

CO 


This condition is almost the same as that in equation (4-3), certainly of the same order of 
magnitude. Thus, collisions are avoided at whole integral periods by the same method as 
they are avoided at half-integral periods; namely, by making the angle <p equal to or 
larger than 0.5 degree. 

If it is desired, at whole integral periods, to keep x larger than some minimum 
value, equation (4-5) yields 


x(nT) = -3(AV cos 0 sin *>)•—- | D • ~~ > x, 

oj 2 co z 


mm 


For horizontal ejections, 0 = 0 degree and cos 0 = 1; thus, the condition on ip is 


cox min , uttD 
-AV sin ip > — + 

6ll7T oo 


or 


< p > sin' 


-i ( ° j2 X m in + 6r2 7f2 D 
l 6n7rcuAV 


(4-7) 


Using typical values of go = 10' 3 , n = 1, D = 10' 6 , and AV = 1 m/s and choosing 
x m ; n = 100 m gives 


ip > sin' 1 (0.01) = 0.5° 


The same considerations for a minimum x-value at one-half an orbital period would 

yield 


p > sin 


, 2aj 2 x min + 3n 2 7r 2 D - 16D 

-1 f 111 111 

\ 6mrcoAV 


(4-8) 


Using approximately the same values gives the same order of magnitude for ip,' i.e., 
approximately 0.5 degree. 


The results, so far, can be summarized as follows: 

1. To keep the TAD’s reasonably close to the Shuttle, within 1 km for example, 
for at least one orbital period, they must be ejected out the side approximately perpen- 
dicular to the orbital velocity vector and within approximately 15 degrees of the local 
horizontal. 

2. To avoid recontact, or collision, exactly sidewise ejections = 0 degree or 
180 degrees must be avoided. To keep the TAD’s at least 100 m away from the Shuttle 
on subsequent passes, the value of should not be less than 0.5 degree. 

Thus, for separation velocities of 1 m/s, the ejection conditions shown in Figure 4 are 
permissible. 

Some examples of relative motions are now presented. Several variations in the 
initial conditions will be made to show how the motions may be varied or how they may 
be controlled through proper selection of the initial conditions. Alt of the motions are 
shown in the relative coordinate system with origin at the Shuttle: +y is along the local 
vertical, +x is opposite the velocity vector, and +z is along the orbital angular momentum 
vector. The motion is three-dimensional. Some two-dimensional cross sections of the 
motion and a three-dimensional plot are presented. 

As a first case four different motions are shown, those of TAD’s ejected into each 
of the quadrants shown in Figure 4. The initial conditions for the four cases are: 


'ase 

0 

JP_ 

AV (m/s) 

a 

+ 5 ° 

+3.5° 

1 

b 

+ 5 ° 

-3.5° 

1 

c 

-5° 

+3.5° 

1 

d 

-5° 

-3.5° 

1 


Figure 5 shows the x-y projection of the motion for these four cases. The solid portions 
of the trajectory correspond to z > 0, and the dotted portions correspond to z < 0. 
Cases b and d, which are forward ejections (negative ip), have motions which carry the 
TAD behind the Shuttle, while cases a and c go in front of the Shuttle. This confirms 
the previous statements. The motions are shown for slightly more than one orbital per- 
iod. The maximum y-displacement during this time is approximately 200 m and occurs 
at approximately one-half an orbital period. At the end of one orbital period, y is back 
to its initial value of zero, x has a secular motion which increases to slightly /more than 
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1 km in one orbital period. In another orbital period the x secular motion would double, 
in another it would triple, etc. The y-motion however would remain periodic except for a 
very slight secular decrease caused by the differential drag. At the orbital altitude of this 
example, 400 km, the differential drag is very small, only about 10' 6 m/s 2 . The motion 
shown in this projection is what an observer traveling parallel with, but at some distance 
off to the side of, the Shuttle would see as the TAD’s were ejected from the Shuttle in 
Iris direction. 

Figure 6 shows the x-z projection of the relative motion for the same four cases. 
This motion is what an observer traveling parallel with, but above, the Shuttle would see 
as he looked down at the Shuttle as the TAD’s were ejected. The solid portions of the 
trajectories correspond to y > 0, and the dotted portions correspond to y < 0. Again, 
cases b and d are observed to go behind the Shuttle, while cases a and c go in front of the 
Shuttle. Case d, which is ejected forward and down ( 6 = -5 degrees, <p = - 3.5 degrees), 
at first goes down below the Shuttle (y < 0) and in front of it. At one-quarter of an 
orbital period it is about 100 m in front of the Shuttle and approximately 900 m off to 
the left side (of an observer looking in the direction of orbital travel). At this point it 
reverses its x-direction of travel and starts in the direction opposite the Shuttle’s direction 
of travel. At the same time y has become positive and comes up above the horizontal 
plane of the Shuttle’s travel. At one-half an orbital period it crosses behind the Shuttle 
at a distance of about 200 m and, as can be seen from Figure 5, about 200 m above the 
Shuttle. This TAD continues traveling in the positive x-direction (opposite the Shuttle’s 
direction of travel). At three-fourths of an orbital period, it is behind the Shuttle about 
800 m and approximately 900 in off to the right side. At the same time its height above 
the horizontal plane is about 200 m, as can be verified from Figure 5. At that point its 
z-motion reverses and it again crosses behind the Shuttle in one orbital period at a dis- 
tance of about 950 m and at the same altitude as the Shuttle (y = 0), as verified by 
Figure 5. If this plot was continued for another orbital period, the motion would repeat 
that of the first orbital period except it would be translated behind the Shuttle about 1 
km. 


The details of the other cases can be determined by inspection of the figures. The 
maximum z-displacement always occurs at one-quarter of an orbital period, the minimum 
z-displacement at three-quarters of an orbital period, and the Shuttle orbital plane (z = 0) 
is crossed at half and whole periods. The oscillatory period of the TAD back and forth 
aecross the Shuttle’s orbit plane is equal to the orbital period of the Shuttle (approxi- 
mately 90 min). There is no way to alter this. 

Figure 7 shows the y-z projection of the relative motion. This is what an observer 
traveling at some distance behind the Shuttle would see as the TAD’s were ejected from 
the Shuttle. The solid portions of the trajectories correspond to x > 0, and the dotted 
portions correspond to x < 0. The symmetric oscillations back and forth across the 
orbit plane are obvious from Figure 7, from which it can be seen that the oscillations in 
y slightly exceed 200 m while those in z are nearly 900 m. The magnitude of these 
oscillations can be controlled (increased or decreased) by changing the magnitude of the 
ejection velocity. 
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The general three-dimensional motion is that of a spiral, although not a very 
symmetric one. The motion generally spirals away from the Shuttle along the x-axis in 
either the positive or negative direction, depending on the ejection conditions. The three- 
dimensional plots make this quite clear. The plots of cases a, b, c, and d above are 
presented in Figures 8 through 11, respectively. The letter code on these plots indicates 
in which octant of the relative coordinate system the TAD is located. The dotted line is 
the projection of the three-dimensional motion onto the x-z plane, i.e., the local hori- 
zontal plane. 

To illustrate the effect produced by changing the ejection velocity, Figures 12, 

13, and 14 are presented. These correspond to Figures 5, 6, and 7 except that the ejec- 
tion velocity has been reduced to one-tenth its initial value, i.e., from 1 m/s to 0.1 m/s. 

By comparing the figures, it can be seen that the motions in each case are similar except 
that the magnitudes have been reduced by approximately a factor of ten. The maximum 
x-displacement in one orbital period is of the order of 100 m rather than 1 km, the maxi- 
mum y-value is of the the order of 20 m rather than 200 m, and the maximum z- 
displacement is about 90 m rather than 900 m. 

To illustrate the effect produced by changing the ejection direction Figures 15, 

16, and 17 are presented. These correspond to Figures 5, 6, and 7 except that the angle 
ip lias been reduced in magnitude from 3.5 degrees to 1 degree. By comparing the x-y 
projections, it can be seen that the maximum y-displacement has been reduced from 200 
m to 100 m and the maximum x-displacement has been reduced from approximately 
1 km to approximately 500 m. The maximum z-displacement is approximately the 
same, however, at approximately 900 m. Comparing the x-z projections, it can be seen 
that the net effect produced by reducing i p is to produce tighter spirals in the x-direction. 

Figures 18, 19, and 20 correspond to Figures 5, 6, and 7 except that the angle 0 
has been reduced in magnitude from 5 degrees to 1 degree. The maximum values of x, 
y, and z are about the same as for the first case. The net result is that cases b and d 
are moved closer together in space and cases a and c are moved closer together. 

To provide more insight into the relative motions, several more variations in the 
initial conditions are shown. Figures 21, 22, and 23 show the x-y, x-z, and y-z projec- 
tions for AV = 1 m/s and (0, <p) = (+10°, +1°), (+10°, -1°), (-10°, +1°), (-10°, -1°). 

These conditions open the spirals considerably, as can be seen, for example, by comparing 
Figure 22 to Figure 16. < p is the same in both cases, but 0 has been doubled in the 
latter instance. Figures 24, 25, and 26 present similiar cases except that 0 has been 
increased even further to 15 degrees and the spirals continue to open further. 

Figures 27, 28, and 29 show cases where y? is held constant at -1 degree and 0 
is varied at 1, 5, 10, and 15 degrees. This distributes the TAD’s at varying distances 
behind the Shuttle. The four cases cross behind the Shuttle at one-half an orbital period 
at distances of 200 m, 450 m, 750 m, and 1050 m, respectively. At one orbital period 
they all sweep behind tire- Shuttle at the common distance of approximately 250 m ' 
providing a good opportunity for contact between two or more of the TAD’s. This con- 
tact probably would not cause any damage, however, since they all have common veloc- 
ities such that relative velocities between them would be near zero. 
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Figures 30, 31, and 32 show cases where 0 is held constant at +5 degrees and 
is varied at -3, -1, +1, and +3 degrees. This distributes the TAD’s both in front of and 
behind the Shuttle, although not in a uniform manner. At one-half an orbital period 
the first three cases cross behind the Shuttle at distances of 750 m, 450 m, and 150 on 
respectively, and the fourth case crosses in front of the Shuttle at a distance of slightly 
more than 100 m. At one orbital period the first two cases cross behind the Shuttle at 
distances of 800 m and 250 m, respectively, while the last two cases cross in front of 
the Shuttle at distances of 350 m and 900 m, respectively. 

Some interest has been expressed in trying to keep the TAD’s exactly in the 
horizontal plane, i.e., no y-component to the motion. Figures 33, 34, and 35 show cases 
where the TAD’s are ejected in the horizontal plane 0 = 0 degree with <. p varied at -3, 

-1, +1, and +3 degrees. Figure 33 shows that although the TAD’s were ejected hori- 
zontally, the subsequent motion is not constrained to the horizontal plane. It is not 
possible to constrain the motion to the horizontal plane because of the coupling of the 
x- and y-motions, i.e., any x-motion will induce a y-motion. 

It is also of interest to determine the effect of increasing atmospheric drag on the 
relative motions of the TAD’s. Figures 36, 37, and 38 show this effect, where the altitude 
has been reduced from 400 km to 300 km and the differential drag has increased by an 
order of magnitude from 1 X JO' 6 m/s 2 to 1 X 10' 5 m/s 2 . These figures are similar to 
Figures 12, 13, and 14 except for the change in drag. The effect produced by the 
increased drag is very pronounced in all the figures. All cases, regardless of variations in 
the initial conditions, show similiar behavior. They quickly decay below the Shuttle and 
rapidly move out in front of it. At one orbital period after ejection, all cases are more 
than 400 m in front of the Shuttle and rapidly moving away. It is doubtful that TAD’s 
would be of much use as a diagnostic instrument for differential drags this high. If it is 
necessary to use them at this altitude (300 km), they probably could be designed small 
and heavy to reduce the differential drag to more acceptable levels. 

Figures 39, 40, and 41 show the same cases at an altitude of 200 km where the 
differential drag has increased another order of magnitude to 1 X 10 -4 m/s 2 . It appears 
that the TAD’s would be completely unuseable at this altitude because of the enormity of 
the differential drag. The TAD’s get 5 km in front of the Shuttle and decay by about 
1 km in altitude in just one orbital period. 

All of the motions shown so far have been for only one orbital period. Figures 
42, 43, and 44 show one case (AV = 1 in/s, 0 = +5 degrees, <p = 1 degree at 400 km 
altitude) for two orbital periods. These figures show the repetitive nature of the motion 
with a slow secular increase in the x-direction. For example at one-half and one and one- 
half orbital periods, the TAD crosses behind the Shuttle at 450 m and 650 m, respec- 
tively, while at one and two orbital periods it crosses at 250 m and 400 m, respectively. 

By proper selection of the initial conditions, it appears that, with a number of 
TAD’s, the region around the Shuttle could -be examined closely. The choice of the ' 
initial conditions would depend on the specific application that one had in mind. 
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V. "MIRROR-IMAGE” MOTION 


If a TAD is ejected with a certain AV, 0, and <p, it has a certain motion relative 
to the Shuttle. If the TAD should then be ejected in the opposite direction (i.e., with 
the same AV but with 0 = -0 and <p = <p +7 r), the resultant motion would be the above 
motion reflected through the origin, i.e., x -> -x, y -> -y, and z -> -z. This is more properly 
termed a coordinate inversion. This can be demonstrated from the equations of motion 
(2-33) and (2-40) with the initial conditions of Section III, 


x = -(AV cos 0 sin ip) (3t - --- sin cot V — Dt 2 + 2 { ~ )(1 - cos cot) 

co / 2 \ co co 2 


y = - “- (AV cos 0 sin ip)(l - cos cot) - — - + f AY, S 2JL?- + ~ } sin cot 


co 


co 


co co 


z = - (AV cos 0 cos i p) sin cot 
co 


(5-1) 


Now, if the replacements 0-^-0 and ->■ <p + n are made, 


sin Op + it) = -sin <p , cos (ip + 7 r) = -cos <p 


sin (- 0 ) = -sin 0 , cos (- 0 ) = +cos 0 


If the terms in equations (5-1) containing D are ignored (since they are small for 
short time periods) and 0 and <p are replaced by -0 and ip + it, respectively, then it 
is easily verified that 


x(-d, tp + 7T) = -x(0 , <p) 
y(-0 )( p + 7r) = -y(d, ip) 


z (-0, <p + 7T) = -z (0, < p) 


(5-2) 


The implication of this is, of course, that ejections on only one side- of the Shuttle, for ' 
example the +z side, have to be analyzed. The motions of TAD’S ejected on the other 
side would be merely a reflection of one of these motions. 
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VI. SENSITIVITY COEFFICIENTS - FIRST ORDER PARTIALS 


Frequently it is of interest to know how sensitive the motion of a TAD is to its 
initial conditions. A practical question arises: Suppose the initial conditions are in error 
by a certain amount because of inaccurate hardware alignment or some other cause, how 
much then is the resultant motion of the TAD affected? Fortunately, this question can 
be answered rather simply for the relative motion equations. The expressions for the 
coordinates are analytic functions of the initial conditions and the time and, hence, can 
be expanded in Taylor series about the initial conditions; i.e., 


x(0,<p,AV,t) 


x(0 o ,ip o ,AV o ,t) + 


(IjA (0-0 o) + 
\o(j J o 



bp - Po ) 


+ 



(AV - AV 0 ) + . . . 


or 


x = x 0 



A(AV) + . . . 


( 6 - 1 ) 


and similiarly for the other coordinates. 

If it is now assumed that the displacements A0, At p, and A(AV) from the nominal 
are small, it is justifiable to retain only first-order terms in the Taylor expansions. The 
resultant variations in the coordinates at any time point resulting from variations in the 
initial conditions will be 



t* 
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These partials are dependent on the nominal values of the initial conditions and, most 
significantly, they are dependent on the time. At some points along the trajectory the 
error in the coordinates because of errors in the ejection conditions will be larger than at 
other points because of time dependence. 

The percentage error in z at least can be shown to be independent of time, 
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Az 

z 


= -tan 0 AO - tan y? Ay? + 


A(AV) 

AV 


(6-4) 


For side ejections where 0 and y? are both small, the percentage error in z is almost 
equal to the last term, the percentage error in AV . 

Some plots of these partial derivatives are presented to get an idea of the numerical 
values they have. Figure 45 shows the partials of the coordinates with respect to 0 for 
one orbital period for the nominal case: 0 = +5 degrees, y? = +3.5 degrees, and AV = 0.1 
m/s. Both 0 and y? are small for this case so the partials of x and y behave like 
the second terms in equations (6-3). i.e., the 3 x/90 dependence on time behaves like 
(1 - cos cot) and that for 9y/9 0 behaves like sin cot. The z-partial remains very small. 
The partials of the coordinates with respect to y? are presented in Figure 46 for one 
orbital period for the same nominal case and the partials of the coordinates with respect 
to the initial impulse AV are presented in Figure 47. 

The following conclusions may be drawn from these figures: 

1. The coordinates are more sensitive to errors in y? than those in 0 by a 
factor of 4 or 5. 

2. The coordinates are very sensitive to an error in the initial impulse, AV. 

3. The largest errors in z occur at one-fourth and three-fourths of an orbital 
period. No error occurs at half and whole orbital periods. 

4. The largest error in y occurs at one-half orbital periods. No error occurs at 
whole orbital periods. 

5. As evidenced by the graphs and the equations, the error in x grows 
monotonically with time. The errors in y and z are periodic in time. 

Figure 48 presents deviations in the coordinates because of “errors” in the ejec- 
tion conditions for the nominal case just discussed. The “errors” are A 0 = 0.5, Ay? = 

0.5 degree, "and A(AV)'= 0.1 m/s. It can be observed that these errors in the coordinates 
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behave exactly like the partials of the coordinates with respect to AV that were pre- 
sented in Figure 47, This is because the coordinates are much more sensitive to errors 
in AV than to errors in 0 or <p, at least for the nominal case being considered here and 
for any nominal case where 0 and <p are both small as they are likely to be for many 
TAD’s applications. 


VII. TARGETING 


There may be applications of the TAD’s for which it would be desirable to place 
them at a certain point in space relative to the Shuttle at a specified time. A specific 
example would be the case where a balloon had been ejected from the Shuttle and a 
TAD would then be ejected such that the TAD would pass through the wake of the 
balloon. We refer to this problem as the targeting problem and formulate it more pre- 
cisely as follows. The TAD is initially in the Shuttle (coordinates: 0, 0, 0,). It is desired 
to have the TAD pass through the point x, y, z in a certain time interval t after ejec- 
tion. What initial velocity x 0 , y 0 , z 0 (or AV, 0, i p) is required to meet these criteria? 

The relative motion equations (2-33) and (2-40) can be written as: 

x(t) = x 0 + 6coty 0 - - Dt 2 - 6y 0 sin cot + (1 - cos cot) + ( -3t + — sin cot) x 0 

2 co^ co 

o 

+ — (1 - cos cot) y 0 

CO 


y(t) = 4y 0 - 



t- 3y 0 cos cot + 



2 

sin cot - (1 

co 


. , . , sin cot . 

- cos cot) x 0 + — — - y 0 
co 


z(t) = z 0 cos cot + — sin cot . (7-1) 

co 

These equations can be rewritten collecting coefficients of x 0 , y 0 , and z 0 on the left 
side with all other terms on the right side. The results are: 


- (-3cot + 4 sin cot) x 0 + (1 - cos cot) y 0 = x(t) - x 0 - 6y 0 (cot - sin cot) 

co co 


+ ~ Dt 2 + ^'5 (1 - cos cot) 


CO 


(7-2) 


- (1 - cos cot) x 0 + - sin cot yo - y(0 - 4y 0 + 2D - t + 3y 0 cos cot - — sin cot 
co co co co 2 
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— sin cot z 0 = z (0 • z o cos cot 
co 


This can now be cast in matrix form as 


1 9 

— (-3cot + 4 sin cot) — (1 - cos cot) 0 

co co 


co 


( 1 - cos cot) 


1 • 

— sm co t 
co 


\ 


Yo 


— sin cot / \ z 0 

CO I 


/ 


3D1 2 4D 

x(t) - x 0 - 6y 0 (cot - sin cot) + + --- (1 - cos cot) 

2 co“ 


?r> on 

y(t) - 4y 0 + — t + 3y 0 cos cot - ~ sin cot 

CO CO 2 


\z(t) - z 0 


cos cot 


(7-3) 


Equation (7-3) is compactly represented by 


Ax 0 = b i (7-4) 

where A is the 3X3 square matrix explicity presented in equation (7-3) and x 0 and 

b are the column matrices in equation (7-3). One inserts the initial position x 0 , y 0 , z 0 
and the desired position at time t, i.e., x(t), y(t), z(t). Equation (7-4) is then solved to 
determine what initial velocity is required to get from one point to another in the speci- 
fied time interval. The solution to equation (7-4) is, of course, given by 


x 0 = A' 1 b 


(7-5) 


The inverse of A in this simple case can be computed by the co-factor method. The 
determinant of A, required by this method, is 


det A = ~ - [-3cot sin cot + 8 (1 - cos cot)] 

or • ■ • 


(7-6) 
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A" 1 exists only if or when det A =A 0. It is easily determined from equation (7-6) that 
the determinant is zero at half and whole orbital periods, so the inverse does not exist 
there. 


A basic reason behind, and the meaning of, the singularity at half and. whole 
orbital periods is that z will be zero at these times regardless of what is done to the 
initial conditions, and it cannot be forced to be otherwise as long as we are restricted to 
a single initial impulse. 

There are also other times at which the determinant is zero, namely those which 
satisfy the transcendental equation, 


f(t) = -Scot sin cut + 8(1 - cos cot) = 0 


(7-7) 


A graph of this function for a few orbital periods will give an estimate of the positions 
of its zeroes. This is shown in Figure 49. Zeroes of this function occur at whole orbital 
periods. In addition, after the first orbital period, there is another zero that occurs just 
prior to a half-orbital period. As t increases this zero more closely approaches one- 
half orbital period. For large t it approaches the zeroes of sin cot , i.e., at half and 
whole orbital periods. This is because the first term in equation (7-7) soon dominates 
the function. 

Thus, det A is zero at half and whole orbital periods and, for the first few orbits 
after the initial one, it is zero at an isolated point just prior to one-half orbit. To calculate 
the position of this zero accurately, equation (7-7) would have to be solved iteratively. 

In practical applications of the targeting procedure, this will cause no difficulty because 
it is extremely unlikely that one would choose a time to target at precisely the position 
of this zero. In summary, then, one can target to a point at any time except at half and 
whole orbital periods. 

With these limitations, the inverse of A is: 


A ' 1 =- 


-3cut sin wt + 8(1 - cos cut) 


j sin cut -2(1- cos cut) 

2(1 - cos cu t) (-3cu t + 4 sin cut) 

^ 0 0 


0 ) 

0 

-3cut sin cut + 8(1 - cos cu t) 
sin cu t j 


(7-8) 
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Multiplying A' 1 by the column matrix b will give the required initial velocity x 0 , y 0 , 
and z 0 . The results of performing this multiplication are: 


x 0 = 


-3wt sin cot + 8(1 - cos cot) 


sin co t 


3 4D 

x(t) - x 0 - 6y 0 (cot - sin cot) + -Dt 2 + — (1 - cos cot) 


2D °D 

-2(1 - cos cot) I y(t) -4y 0 + — t + 3y 0 cos cot - ~ sin cot 


Yo = 


■3co t sin co t + 8( I - cos co t) 


jj2( l - cos cot) | ; 


. 4D. 


2(1 - cos cot) { x(t) - x 0 - 6y 0 (cot - sin cot) + -Dt 2 + — (I - cos cot) 

9 . t. j 


9Q 2D 

+(-3cot + 4 sin cot) { y(t) - 4y 0 + — t + 3y 0 cos cot - ■ sin cot 

to co 1 


z 0 - . -~ w — [z(t)-z 0 coscot] , (7-9) 

sin co t 

It can be seen from these equations that x and y are coupled but z is not. If z is 
required to be other than zero at half and whole orbital periods, the required z 0 would 
be infinite, which is an impossibility. If x and y were required to be some arbitrary 
value at one orbital period, the y 0 equation, from (7-9), would be 


6cott 


Vo 


,,, , 4ttD 

yCt) - y 0 + -- 

ojr _ 


(7-10) 


This is infinite unless the term in brackets is zero; i.e., unless 


y(t) 


Yo ' 


47tD 

oj 2 


(7-1 1) 


Thus, the only specification that can be made on y at one orbital period is that it 
returns to its initial value minus the amount that it decays in one orbit because of drag. 

The x 0 equation yields the indeterminant form 0/0 at one orbital period. An 
application of L’HospitaTs rule with subsequent evaluation at one orbital period yields 


gj(x -x 0 ) , n ttD 

- - - - 0 + 2wy 0 - - 

■ -07T ' ■ ' -- GcJ • 


(7-12) 
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Exactly the same result can be obtained from equations (2-40) and (2-41). Since, in our 
case, x 0 and y 0 are both zero, equation (7-12) becomes 


Thus, we get the important result that although y and z cannot be specified at one 
orbital period (they both return to zero), x can be specified at one orbital period and 
equation (7-13) gives the initial impulse required (the x-component). 

It can similiarly be determined that although z cannot be specified at half-orbital 
periods, both x and y can be. The results can be summarized as follows: 

1. One can target to an arbitrary point (x, y, z) at times other than half and 
whole orbital periods. 

2. One can target to an arbitrary (x, y) point (but not z) at half-orbital periods. 

3. One can target to an arbitrary x-point (but not y or z) at whole orbital 

periods. 

Once the target point (x, y, z) and the time t are specified, equations (7-9) determine 
the required velocity. The required AV, 6 , and y will be 


AV = J x 0 2 + Yo 2 + Zo 2 


= sin' 1 / 

\ AV/ 


i / *0 


These equations are evident from Figure 3. 


Two numerical examples of the targeting capability are presented. The two 
examples are: target to the point x = 200 m, y = 200 m, z = 200 m in (a) one-fourth 
orbital period and (b) one and one-fourth orbital periods. The required initial conditions 
from equations (7-14) are: 
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Case a 


x 0 = -0.06844558 m/s 

AV = 0.253099 m/s 

y 0 = +0.0903975 m/s 

0 = 20.926° 

z 0 = +0.22628 m/s 

<p = -16.8296° 

Case b 

x Q = +0.0108422 m/s 

AV = 0.3449086 m/s 

y 0 = +0.26008036 m/s 

0 = 48.9429° 

z 0 = +0.22628 m/s 

p = 2.7432° 


The x-y projection of the motion for both cases is presented in Figure 50. Case a is 
shown for approximately one orbital period and case b is shown for approximately two 
orbital periods, The x-z projection is presented in Figure 51 and the y-z projection is 
presented in Figure 52. It can be verified that both cases go through the target point at 
the specified time. 

In order to hit the specified target point within a certain accuracy, it is necessary 
to determine the accuracy with which the initial velocity must be controlled. This can 
be determined with the aid of the partials developed in Section VI. For case a the partials 
of the coordinates with respect to 0 are presented in Figure 53 for one orbital period. 
The partials of the coordinates with respect to i p are presented in Figure 54 and the 
partials of the coordinates with respect to AV are presented in Figure 55. From these 
figures it can be determined at one-quarter of an orbital period, the time at which we 
desire to hit the target point, that the partials are approximately 
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II: seems remarkable that the partials of the coordinates with respect to the initial AV 
are all the same at one-quarter period. IT the allowable errors Ax, Ay, and Ax in the 
position x, y, z are specified, the allowable errors in the ejection conditions can be 
determined by inverting the matrix in equation (6-2). If equation (6-2) is written as 


Ax(t) = B(t) AVj 


(7-16) 


then 


AV; = B" 1 (t) Ax(t) 


(7-17) 


It can be verified, by much tedious algebra, that the inverse of the partial matrix 
B, defined by 
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where 


|B| = [-■ (- cos 0 sin cot: [Scot sin cot - 8(1 - cos cot)] 

AV \ co / 


is 


Bit = 


AV 

b, , - ~ r 


/av ' 2 

12 AV \co 


J cos 0 sin cot [sin 0 sin sin cot - 2 cos 0(1 - cos cot)] 

cos 0 sin cot (.^^.sin cot 
( cos 0 


(7-19) 


B, 3 = ^ cos 0 sin cot [-2 sin 0(1 - cos cot) - sin y cos 0 sin cot: j 

Bn = ° 0S ^ S ' n ^ s * n co ^ cos 0 - Z sin 0 sin </>(! - cos cot) j 
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The zeroes of the determinant |B| occur at the same places as those for the determinant 
|A|, i.e., at half and whole orbital periods; therefore, B' 1 does not exist there just as 
A" 1 does not. In addition, B" 1 does not exist for vertical ejections, i.e., for 9 = tt/ 2 or 
9 ~ 37r/2 . 

We divide each element in equation (7-18) by |B| and, for example, rewrite 
B U /|B| = bji. Equation (7-17) then becomes 
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where each element can be explicitly determined from equations (7-19) and (7-20). The 
units on each of the elements in the first two rows are radians/meter. These can be 
multiplied by 57.29577951 to change them to the more convenient units of degrees/meter. 
The units on the elements in the last row are (seconds) -1 . The elements in the first row 
of equation (7-21) are plotted in Figure 56 for case a for one orbital period. The units 
are shown as degrees/meter. The elements in the second row are plotted in Figure 57 and 
the elements hi the third row are plotted in Figure 58. 

The elements at one-quarter period, reading from these graphs, are approximately 
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If we suppose that up to a +10 m error could be tolerated in all of the coordinates at 
this time, the allowable errors in the ejection conditions would be 


AO = -0.05° 

A<p = -0.0 6° 

A(AV) = 0.0125 m/s 


If these errors are multiplied by the partial matrix in equation (7-15) it is verified that 
they will indeed produce approximately a +10 m error in all of the coordinates. If up 
to a -10 m error in all coordinates could be tolerated, the errors allowed in the ejection 
conditions would be the negative of those above. If a +10 m error is allowed in x and 
none in y or z, the allowable errors in the ejection conditions would be AO = +1.5°, 

Aip = +0.77°, and A(AV) = +0.0016 m/s. 

One aspect becomes clear from this; that is, if an allowed “box” of errors in the 
coordinates at some point in time is specified, this can be mapped into an allowed volume 
of errors in the ejection conditions. This provides a method by which one may specify 
how accurately the ejection conditions 0, <p, and AV must be controlled. 

The smallness of the coefficients in the third row of equation (7-22), and from 
Figure 58, indicates that the impulsive velocity will have to be controlled very accurately. 
For the example quoted above, it cannot have an error of more than 0.01 m/s. 

As a concrete example of the mapping capability, a cube of allowed position 
errors having sides of 20 m in magnitude centered about the nominal point is taken at 
20 min after ejection and mapped back into allowed ejection errors in 0, tp, and AV. 

The A0, Aip results are plotted in Figure 59. The allowed error in AV, while not 
shown on Figure 59, is always equal to less than 0.01 m/s in absolute magnitude. It is 
seen that the allowed “box” of ejection errors is a rectilinear figure. This is because the 
mapping, being linear, maps straight lines into straight lines. Thus, from this figure, one 
sees that the allowed error in 0 can be up to approximately 3 degrees and that in ip 
up to approximately 3.5 degrees, and the TAD will still be somewhere in the cube having 
sides of 20 m defined about the nominal point at 20 min after ejection (provided that 
the error in AV is also less than 0.01 m/s). At various time points along the trajectory 
this allowed error box will vary some in size and shape. As an example, Figure 60 gives 
the allowed error box at 25 min. Although they look similiar at first glance, a careful 
comparison will reveal differences. 

Another important point should be emphasized about equation (7-21); that is, 
the allowed ejection errors are linear functions of the allowable errors in the final position, 
so the allowable ejection errors scale linearly up or down with the allowable position 
errors. 
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The targeting capability can also be used to maximize the amount of space sampled 
by a given number of TAD’s. For example, suppose it was desired to examine the wake 
behind the Shuttle with four TAD’s. These TAD’s could be distributed, at one-half an 
orbital period or at a whole orbital period, at the best distances as determined by the 
investigators; they could be distributed at 250, 500, 750 and 1000 m or at any other set 
of desired distances. If it was desired to make them pass through the wake, not simul- 
taneously but in a certain time sequence, this could easily be accomplished by proper 
time-sequencing of the ejections. 

The discussion concerning targeting accuracy thus far has assumed that the differ- 
ential drag is known perfectly. In fact there could be, and probably will be, some uncer- 
tainty about its exact value and this could upset the accuracy estimates. To determine 
the effects that an uncertainty in drag would produce, it is necessary to look at the 
partials of the coordinates with respect to the differential drag. From equations (7-1) it 
is easily determined that these partials are 
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The z-coordinate is not affected at all by drag or by any uncertainty in it. It is seen that 
these partials are not dependent on initial conditions and only slightly dependent on 
altitude (through co). Therefore, they are approximately the same for all orbits. A plot 
of these partials versus tune for one orbital period is given in Figure 61. It is seen that 
the secular terms quickly dominate the periodic ones, and the partials decrease monoton- 
ically with time indicating that the errors introduced by uncertainties in drag would be 
larger for large times than for small tunes. This is what one would intuitively expect. 

Ignoring other uncertainties for the present, the errors induced in the coordinates 
because of an uncertainty AD in the differential drag would be 
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The uncertainty AD in differential drag could be constant or time varying. If it was time 
varying, however, its exact form would probably be unknown to us. In the simple case 
where AD is taken to be a constant value, the errors in the coordinates at any time 
point are obtained simply by multiplying the partials at that time point by the uncertainty. 

At an altitude of 400 km, the nominal differential drag is of the order of 10' 6 
m/s 2 . If the uncertainty is 10 percent of the nominal value, AD would be of the order 
of ±10" 7 m/s 2 . In this case the induced errors in the coordinates over one orbital period 
would be negligible (< 1 m). If the uncertainty is 100 percent of the nominal value, the 
induced errors in the coordinates could be ±10 m or more in the first orbital period (and 
larger in succeeding orbital periods). This may be significant for precise requirements. At 
an altitude of 300 km where the nominal differential drag is of the order of 10“ s m/s 2 , 
one would definitely have to become concerned about uncertainties in differential drag. 


VIII. OPTIMUM USE OF TAD's FOR ELECTROMAGNETIC INTERFERENCE 

AND WAKE STUDIES 


One of the desired objectives of this study is to identify trajectories of TAD’s 
that will cover or scan as much volume around the Shuttle as possible. It is desired to 
maximize the volume coverage while minimizing the number of TAD’s used. Since 
TAD’s can be made to go forward and backward and to oscillate in the vertical plane 
and sweep back and forth in the horizontal plane (although these do not remain in the 
horizontal plane but have vertical displacements), it would seem at first glance that four 
TAD’s would give reasonably complete coverage of the volume around the Shuttle. The 
plan is to have two TAD’s oscillating in the vertical plane with one in front of the 
Shuttle and one behind the Shuttle. The other two TAD’s would be given horizontal 
sweeping motion with one in front of the Shuttle and one behind the Shuttle. 

With this plan in mind and assuming that 100 m above and below and to each side 
of the Shuttle is a sufficient distance for electromagnetic interference and wake studies, it 
is determined from equations (4-1) that an impulse of 0.1 m/s is sufficient to achieve 
these distances in oscillations. For the vertical oscillations, the angle i p should be 90 
degrees (in-plane ejections). For the TAD going behind the Shuttle, it can be determined 
that a 0 angle of 100 degrees (slightly forward of vertical) will produce an acceptable 
trajectory and for the TAD going in front of the Shuttle, a 0 angle of 270 degrees 
(straight down) will produce an acceptable trajectory. 

These two trajectories are shown in Figure 62 for approximately five orbital 
revolutions. Case a going in front starts with overlapping trajectories but drag starts 
drawing them out until by the fourth orbit they are no longer overlapping. Drag is 
also slowly making them decay at a mb' of about 10 m per orbit (this is at an altitude 
of 400 km with a differential drag of 1 X HT 6 m/s 2 ). This TAD will continue to acceler- 
ate away from the Orbiter. Case b going behind the Shuttle starts receding from the 
Shuttle, but by the third orbit drag has stopped the receding motion at about 800 m 
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and starts the TAD back toward the Shuttle. In the five orbits shown here, the TAD has 
not yet gotten back to the Shuttle but will in one or two more orbits. This gives a 
definite possibility of recontact and is a potential safety problem. The relative velocities 
between the TAD and Shuttle are so small (approximately 0.1 m/s), however, that even 
if recontact did occur it is doubtful that it would result in any damage. It might even be 
possible to recover the TAD if it came close enough. (The continuation of this trajectory 
is shown later.) These trajectories make it appear that the vertical plane in front of and 
especially that behind the Shuttle can be sampled quite densely by TAD’s for several 
orbital periods. 

For the two TAD’s executing horizontal sweeps, the initial angle 0 is chosen 
to be zero. For the TAD going forward, the initial <p is chosen to be 2 degrees and for 
the TAD going backward, the initial < p is chosen to be -12 degrees. These two cases, 
labeled c and d, are shown in Figure 63 for approximately five orbits. The vertical plane 
(x-y) motion for these two cases is shown in Figure 64. Case c going in front of the 
Shuttle starts off with a compact oscillation, the distance between successive axis crossings 
being about 100 m. Just as in case a previously, however, the differential drag stretches 
these oscillations and this TAD is accelerated away from the Shuttle. The TAD in case 
d going behind the Shuttle is rapidly slowed by drag and its backward motion is halted 
at about 3.75 orbital periods, at which time it starts forward coming back toward the 
Shuttle. It is fortuitous that the backward motion in this case ceases at almost exactly 
3.75 orbital periods, because this makes the motion coming back toward the Shuttle 
follow almost the same projected path in the horizontal plane as it took going out. This 
poses a definite safety problem. (The continuation of this trajectory is also presented in 
the following.) 

The continuation of case b, initially presented in Figure 62, is shown for about 
nine orbits in Figure 65. The TAD passes back by the Shuttle on the sixth orbit at 630 
min after ejection. At this time it is about 130 m below the Shuttle so there is no recon- 
tact problem in this case. After the TAD passes the Shuttle, it is then accelerated away 
and there is no further concern with recontact. The closest approach of this TAD to the 
Shuttle as it returns occurs at about 565 min, where it crosses the x-axis immediately 
behind the Shuttle at a distance of about 40 m. Any slight variation in the conditions 
of this case, such as differential drag or ejection conditions, could critically alter it, 
changing it from a safe trajectory to an unsafe one from a recontact point of view. 

The continuation of case d, initially presented in Figure 63, is shown for about 
ten orbits in Figure 66, with the vertical motion shown in Figure 67. The horizontally 
projected oscillations coming back follow the same path as they went out on and go 
right back through the origin. The vertical projection of the motion coming back is the 
mirror image of that going out. The projection, too, goes through the origin. This trajec- 
tory would definitely be a problem unless a provision was made to “catch” the TAD as 
it came back or to move the Shuttle to a different position before the TAD got back. 

One nice feature about these trajectories is that double coverage of the wake 
behind the Shuttle is obtained, i.e., going out and coming back. As the TAD passes by 
coming back, it might be possible to recover and reuse it. Another possibility is to design 
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the trajectory such that the TAD misses the Shuttle as it comes back and then continues 
on to map the region in front of the Shuttle. With this possibility, the region behind and 
in front of the Shuttle could be mapped with one or two TAD’s rather than four. For 
completeness the y-z projection of this motion is shown in Figure 68. 

To show that the horizontal projection of the returning motion following that of 
the outgoing motion was coincidental in the preceding case, we present a similiar case 
where this does not happen. The ejection conditions are changed slightly. The angle <p 
is changed from -12 degrees to -14 degrees. All other values remain fixed. The motion 
in this case is presented in Figures 69 through 71. The projection of the returning motion 
here clearly does not follow that of the outgoing motion. In this case the returning TAD 
easily misses the Shuttle. This TAD could give double coverage of the wake behind the 
Shuttle to about 1 km and then could map some of the region in front of the Shuttle. 

The question remains: How does one choose these trajectories optimally? To 
answer this quantitatively it would be helpful if some specific measure could be chosen 
to optimize. Then the pertinent parameters can be chosen so as to optimize that measure. 
As a specific example of this possibility, we might choose to maximize the amount of 
time that the TAD spent behind the Shuttle. To do this we take the secular terms in the 
x-equation (4-1) and equate them to zero. (The oscillatory terms will average out to zero 
and hence can be ignored here.) This yields 


, _ -3AV cos 0 sin ip 

* 1 .5 D (8 -° 

as the amount of time spent behind the Shuttle. This can be made as large as desired by 
choosing AV as large as desired. This, however, would not satisfy our desire to remain 
close to the wake of the Shuttle. Obviously we have constraints on the optimum problem; 
we must maximize t while remaining close to the wake of the Shuttle. If we want to 
remain within ±100 m of the wake, AV must be constrained to approximately 0.1 m/s. 

D is generally fixed by the orbital altitude and the configurations and attitudes of the 
spacecraft involved, so the optimization must be through manipulation of the angles 0 
and ip. If we hold 6 constant at zero and vary < p, we get the result shown in Figure 
72. If 0 is varied, plus or minus from zero, the result will be lower than that shown. 
Clearly, then, t is maximized by chosing <p = -90 degrees and 0=0 degree. 

This is one way to optimize the trajectory selection: choose a quantitative meas- 
ure, and maximize or minimize it. The problem is finding the proper quantity to opti- 
mize. There may be undesirable features about maximizing t above. For example, the 
maximum distance that the TAD goes behind the Shuttle in this case is about 15 km 
which may be farther than experimental interests dictate. The x-distanee behind the 
Shuttle versus time is shown in Figure 73 with i p parameterized. From this figure one 
can pick <p such that x does not go beyond a certain distance. 
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One may be interested in having several wake crossings spaced closely behind the 
Shuttle rather than extended out over a large distance behind the Shuttle. Figure 74 
shows the distances of the first five wake crossings behind the Shuttle versus the ejection 
angle i p. For <p = -10 degrees, the first five crossings occur within 0.5 km of the Shuttle. 
For i p ~ -30 degrees the first five crossings are distributed out to 1.8 km, and for ip = 

-90 degrees they are distributed out to 4 km. The local altitude (y-coordinate) above or 
below the wake of the Shuttle at the wake crossings (z = 0) is shown for the first six 
wake crossings versus i p in Figure 75. The altitudes of the odd-numbered wake crossings 
(half-orbital periods) increase rapidly with increasing i p from about 50 m at <p = -10 
degrees to about 350 m at < p = -90 degrees. This is lowered about 10 m per orbit for 
each successive odd-numbered wake crossing because of the relative decay of the TAD 
orbit caused by the differential drag. 

The even-numbered wake crossings (whole orbital periods) would all be zero if 
there was no differential drag. The drag considered here (10" c m/s 2 ) causes each orbit 
to decay about 10 m below that of the previous orbit. 

Information equivalent to that shown in Figure 74 is presented in a different form 
hi Figure 76. This shows the distance of the wake crossings behind the Shuttle versus 
wake crossing number with cp parameterized from -10 degrees to -90 degrees. 

Figure 77 presents the distance of the TAD in front of the Shuttle versus time 
with i p parameterized. This figure complements Figure 73. From this figure one can 
determine, for example, that the TAD ejected at an angle of ip = -10 degrees will stay 
hi front of the Shuttle and within 1 km of the Shuttle for about 4 hours. The TAD 
ejected at ip = -90 degrees will stay within this region for only 1 hour. This might be a 
useful criterion for choosing optimum trajectories; i.e., maximize the amount of time that 
the TAD stays within a certain distance of the Shuttle. 

To illustrate the above possibility, the amount of tune that the TAD is behind 
the Shuttle and within a certain distance of the Shuttle versus the ejection angle i p is 
presented in Figure 78. For example, if we wanted to maximize the amount of time 
that the TAD is within 2 km of the Shuttle, we would go along the line labeled x < 2 
km until we found its maximum point. In this case it occurs at <p = -20 degrees and 
gives 19.5 hours of time (or about 13 orbits) behind the Shuttle with x < 2 km. Any 
other ejection angle will give less time in this region. The height at the first wake crossing 
in this case, from Figure 75, would be about 116 m and that at the second wake crossing 
would be about -10 m. The wake crossings on each successive orbit would be about 10 m 
less; so, by the 12th orbit the altitude of the wake crossing on the half-orbit would be 
slightly negative (-1.7 m from calculations) and the altitude of the wake crossing on the 
whole orbit would be -1 18 m. Thus, this orbit starts with nearly all of the vertical dis- 
placement above the horizontal plane and ends with all of it below the horizontal plane. 
This TAD would pass beneath the Shuttle and continue to accelerate down and away 
from the Shuttle. 


43 



In Figure 79 we present the amount of time that the TAD is within a given dis- 
tance in front of the Shuttle versus the ejection angle <p. These curves all approach a 
maximum for the smallest possible <p. 

Figure 80 gives the amount of time that the TAD is within plus or minus a given 
distance of the Shuttle versus the ejection angle <p. It is a combination of the two pre- 
vious figures. The maximums on it occur at the same ip’s as the maximums on 
Figure 78. 

From the above discussion it can be seen that the following criteria are among 
the probable numerous ones available to pick from to optimize the TAD trajectories: 

1. Maximize the amount of time behind the Shuttle. 

2. Maximize the amount of time within a given distance behind the Shuttle. 

3. Maximize the amount of time within plus or minus a given distance of the 

Shuttle. 

4. Do any of the above while limiting the vertical and horizontal oscillations to 
some given value. 

From the discussion in this section it is known that if there are no constraints on 
maximizing the amount of time behind the Shuttle, that it is practically unlimited. The 
figures that have been presented enable one to pick optimum values for the remaining 
possibilities for the single case where AV has been limited to 0.1 m/s. 


IX. SUMMARY 


Hopefully, it has been shown from an operational point of view that TAD’s can 
be useful as diagnostic tools in the AMPS program. Their motions can be controlled to 
a very significant degree simply through careful control of their ejection conditions. It 
has been shown that they can be made to go through a specified point at a specified 
time, and a technique has been developed to determine how accurately the initial condi- 
tions must be controlled. Those initial conditions that might result in recontact with the 
Shuttle have been identified and those which cause the TAD to rapidly leave the vicinity 
of the Shuttle have been identified. This results in a range of acceptable initial condi- 
tions being delineated. 

The effect of atmospheric drag on the motion of the TAD’s has been demonstrated. 
It appears doubtful that they can be useful at altitudes as low as 200 km because of the 
very strong effect of drag there. They probably can be useful at 300 km and most cer- 
tainly can be useful at 400 km and above. 


The targeting capability may prove to be extremely useful because it will allow 
one to examine whatever point may be of interest to him. The only limitation to examin- 
ing a large number of points would be the number of TAD’s available because, in general, 
it will require a separate TAD for each discrete point to be examined (unless the TAD’s 
lire recovered and reused). 

A brief look was taken at what an uncertainty in differential drag might do to 
targeting accuracy. It was shown that at high altitudes (> 400 km), it is probably of 
little consequence. At lower altitudes (< 300 km), however, it could definitely pose a 
significant problem. 
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Figure 6, The x-z projection of the relative motion. 
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Figure 7. The y-z projection of the relative motion. 
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Figure 8. Three-dimensional plot of the motion of a TAD relative to the Space Shuttle 

for approximately one orbit. 
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Figure 9. T teee-dimensio,i^ plot of the motion of a TAD relative to the 
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400 km CIRCULAR ORBIT D = 1X 10" 6 m/s 2 

INITIAL CONDITIONS: 

CASE c 0 = -5° 

(p = +3.5° 

AV = 1 m/s 

Figure 10. Three-dimensional plot of the motion of a TAD relative to the 
Space Shuttle for approximately one orbit. 


53 


LETTER CODE 
V 



400 km CIRCULAR ORBIT D = 1 X 10' 6 m/s 2 

INITIAL CONDITIONS: 

CASE d 0 = -5° 

<p = -3.5° 
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Figure 11. Three-dimensional plot of the motion of a TAD relative to the 
Space Shuttle for approximately one orbit. 
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Figure 12. The x-v projection of the relative motion. 
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Figure 13. The x-z projection of the relative motion. 


56 



d 


d 


+Z(m) 


i 

i ■ • • \ > 

a 

+100 

c 



b 

c - 100 
a 


-100 


x > 0 
x < 0 


400 km CIRCULAR ORBIT 


D = 1 X 10' 6 m/s 

CASE 

0 (degrees) 

i p (degrees) 

AV (m/s) 

a 

+5 

+3.5 

0.1 

b 

+5 

-3.5 

0.1 

c 

-5 

+3.5 

0.1 

d 

-5 

-3.5 

0.1 


DOTS ARE SEPARATED BY 1 min IN TIME. ORBITAL VELOCITY 
VECTOR INTO PAGE. 

Figure 14. The y-z projection of the relative motion. 
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Figure 16. The x-z projection of the relative motion. 
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Figure 18. The x-y projection of the relative motion 
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Figure 19. The x-z projection of the relative motion 
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Figure 20. The y-z projection of the relative motion. 
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Figure 21. The x-y projection of the relative motion. 
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Figure 22. The x-z projection of the relative motion. 
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Figure 23. The y-z projection of the relative motion. 
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Figure 25. The x-z projection of the relative motion. 
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Figure 26. The y-z projection of the relative motion. 
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Figure 27. The x-y projection of the relative motion. 
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Figure 28. The x-z projection of the relative motion. 
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Figure 29. The y-z projection of the relative motion. 
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Figure 30. The x-y projection of the relative motion. 
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Figure 31. The x-z projection of the relative motion. 
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Figure 32. The y-z projection of the relative motion. 
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Figure 33. The x-y projection of the relative motion. 
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Figure 34. The x-z projection of the relative motion. 
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Figure 35. The y-z projection of the relative motion. 
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Figure 36. The x-y projection of the relative motion. 
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Figure 37. The x-z projection of the relative motion. 
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Figure 38. The y-z projection of the relative motion. 
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Figure 40. The x-z projection of the relative motion. 
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Figure 41. The y-z projection of the relative motion. 
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Figure 42. The x-y projection of the relative motion. 
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Figure 43. The x-z projection of the relative motion. 
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Figure 44. The y-z projection of the relative motion. 
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Figure 45. Partials of the coordinates with respect to the initial angle 0 
versus time for the nominal case ( 0 = +5 degrees, = +3.5 degrees, 

AV = O.i m/s). 
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Figure 47. Partials of the coordinates with respect to the initial impulse 
AV versus time for the nominal case ( d = +5 degrees, 

V? = +3.5 degrees, AV = 0.1 m/s). 
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Figure 49. A plot of the function f(t) = -3cut sin cut + 8(1 - cos cut) 
versus time for three orbital periods. 
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Figure 50. The x-y projection of the relative motion, (two sample cases 
targeting to the point in (a) one-fourth T and (b) one and 
one-fourth T] . 
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Figure 51. The x-z projection of the 
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one-foui 
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Figure 52. The y-z projection of the relative motion [two sample cases 
targeting to the point in (a) one-fourth T and (b) one and 
one-fourth T] . 
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Figure 53. The partials of the coordinates with respect to the initial angle 
9 versus time for the nominal case ( 9 = 20.926 degrees, = -16.829 
degrees, AV = 0.253 m/s). 
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Figure 57. Coefficients determining allowable error in the ejection angle <p 
for a specified allowable error in the position versus time for the 
nominal case (0 = 20.926 degrees, tp = -16.829 degrees, 

AV - 0.253 m/s). 
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Figure 58. Coefficients determining allowable error in the ejection impulse 
AV for a specified allowable error in the position versus time for the nominal 
case (0 = 20.926 degrees , \p - -16.829 degrees, AV = 0.253 m/s). 
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Figure 59. The allowed ejection error box for a cube of allowed position 
errors of sides 20 m about the nominal point at 20 mm for the nominal 
case (9 = 20.926 degrees, ip = -16.830 degrees, AV — 0.253 m/s). 
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Figure 60. The allowed ejection error box for a cube of allowed position 
errors of sides 20 m about the nominal point at 25 min for the 
nominal case ( d = 20.926 degrees, = -16.830 degrees, 

AV = 0.253 m/s). 
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Figure 62. The x-y projection of the relative motion. 
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Figure 63. The x-z projection of the relative motion. 








400 km CIRCULAR ORBIT D = 1 X 10‘ 6 m/* 2 

CASE b: Q = 100°, » 90°, AV = 0.1 m/s 

THIS CASE SHOWN FOR APPROXIMATELY NINE ORBITS 
DOTS ARE SEPARATED BY 10 min IN TIME 
ALL MOTION IN THE ORBITAL PLANE (x-y) FOR THIS CASE 

Figure 65. The x-y projection of the relative motion. 
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Figure 66. The x-z projection of the relative motion. 
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Figure 67. The x-y projection of the relative motion. 
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Figure 68. The y-z projection of the relative motion. 
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Figure 69. The x-z projection of the relative motion. 


2 


0 km CIRCULAR ORBIT 


D = 1 X 10*° m/* 


CASE e: 0 = 0°, <p = -14°, AV = 0.1 m/s 

■ . 

THIS CASE SHOWN FOR APPROXIMATELY TEN ORBITS 
DOTS ARE SEPARATED BY 10 min IN TIME 

Figure 70. The x-y projection of the relative motion. 





Figure 72. Time TAD stays behind the Shuttle versus the 
ejection angle i p. 
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Figure 74. Distances of the first five wake crossings behind the Shuttle 
versus the ejection angle i p. 
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Figure 77. Distance of the TAD in front of the Shuttle versus time with 
the ejection angle \p parameterized. 
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Figure 78. Time of TAD behind the Shuttle within a given distance versus the 

ejection angle ip. 
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Figure 79. Time of TAD in front of the Shuttle within a given distance 

versus the ejection angle <p. 
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Figure 80. Time of TAD within plus or minus a given distance of the Shuttle 

versus the ejection angle <p. 
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